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The lowest-energy structures of lead clusters interacting via a Gupta potential are obtained for 
TV < 150. Structures based on Marks decahedra dominate at the larger sizes. These results are very 
different from those obtained previously using a lead glue potential, and the origins of the differences 
are related back to differences in the potential. 
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I. INTRODUCTION 



II. METHODS 



The structure of a metal cluster is one of its most im- 
portant properties, yet it is a property that can be hard 
to access for free clusters. Experiments usually lead only 
to an indirect measurement of the structure, and so re- 
sults need to be compared to what would be expected 
for candidate structures. Therefore, there is an impor- 
tant need for reliable predictions from theory. However, 
this is also something that can be hard to achieve. 



For all but quite small clusters, performing a full search 
for the lowest-energy structure is unfeasible with ab ini- 
tio electronic structure methods. Therefore, one needs 
to use some kind of potential. However, the energy dif- 
ferences between competing structures can often be rel- 
atively small. Furthermore, for a cluster's energy to be 
accurate, the energetics associated with a wide range of 
properties potentially need to be well modelled — different 
external surfaces, twin planes, different crystal struc- 
tures, and the response to strain. To model the poten- 
tial temperature dependence of the structure vibrational 
properties need also to be well described. 1 Therefore, pre- 
diction of the correct structure of a cluster represents a 
tough challenge for a potential. Furthermore, often it is 
not clear which of the available potentials for a system 
will work best. 



It is my contention that an important aspect of as- 
sessing a potential's quality for modelling clusters is to 
first understand the relationship between the observed 
structures and the potential. As well as providing fun- 
damental insights into cluster structure, this can help to 
determine whether the observed structures reflect some 
robust feature of the potential, or rather whether they 
stem from some deficiencies in the potential. This is the 
approach I undertake here. I first find the lowest-energy 
structures for lead clusters modelled using a Gupta po- 
tential, and then compare them to previous results for 
a lead glue potential^ The origins of the differences are 
explored, highlighting the implications for the relative 
merits of the two potentials. 



Both potentials that I consider are of the embedded- 
atom form, where the potential energy is given by 
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where cf>(r) is a short-ranged pair potential, F{p) is 
a many-body embedding (or glue) function and pi = 
P { r ij)i where p{r) is an "atomic density" function. 
One potential is of the Gupta form: 3 
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whereas for the glue potential there are no assumed 
forms for these functions; instead they arise in the fitting 
processil The parameters of the Gupta potential were 
fitted on the basis of bulk properties^ whereas a much 
wider variety of configurations, including surface proper- 
ties, were used to parameterize the glue potential. 4 

However, the functions in Eq.|21are non-unique. Func- 
tions that give rise to exactly the same energy can be 
constructed by the transformation 
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This transformation redistributes the total energy be- 
tween £ pair and -Eembed- When 
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F'{p) has a minimum at p x tai, where p x t a i is the value 
of p in the equilibrium crystal. This choice is called the 
effective pair format, and has been suggested as the most 
natural way to partition the energy between the pair and 
many-body contributions^ In this format, when p = p x tai 
the pair potential controls the energy change for any 
change of configuration that does not significantly alter 
p. This form also makes comparisons between different 
potentials easier. 
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FIG. 1: (a) 4>(r) and p{r), and (b) F(p) for lead Gupta and glue potentials in the effective pair format. 



Indeed, the components of the two lead potentials are 
compared in Fig. ^ and are quite different. The pair 
potential contributes less for the glue potential. At its 
minimum the energy is only 1.62% of the minimum in the 
embedding function, compared to 3.94% for the Gupta 
form. The p(r) functions are also very different. For 
the Gupta potential it decreases rapidly and uniformly 
because of its exponential form, whereas the glue p(r) has 
a wide plateau around the nearest-neighbour distance. 

The optimal clusters will have the best balance be- 
tween maximizing the pair and embedding energies. As 
is well-understoodf& maximization of the pair energy 
is achieved through a balance between maximizing the 
number of bonds (i.e. spherical and facets with high co- 
ordinate surface atoms, e.g. fee {111}), and minimizing 
the strain energy (the energetic penalty for distances that 
deviate from the minimum of the potential). 

The way the embedding energy is maximized and the 
effects on structure have only begun to be explored more 
recentlyi 2 ^!?!?! 1 ?! 11 ! 12 For i? cm bod to be maximal the pi 
for each atom need to be as close to p xta i as possible. 
An important feature is that an atom with a low (high) 
coordination number can compensate for this by short- 
ening (extending) the bonds. Therefore, the embedding 
energy will favour a shrinking of the surface, which can 
in turn compress the interior of the cluster. This can be 
problematic for structures such as the Mackay icosahe- 
dron, where bond distances for surface atoms are longer 
than interior atoms, and so surface compression can lead 
to unfavourable energies for the atoms at the core. The 
magnitude of these effects will depend on the relative 
depths of the effective pair and embedding functions, as 
the pair potential will tend to resist any changes in struc- 
ture that lead to distances deviating from the preferred 
pair distance. When the embedding energy dominates 
unusual structures that naturally have bond distances 
that are shorter on the surface can be favourediii*^ 

For Gupta potentials the value of p/q provides an ap- 
proximate indicator of the relative weights of the pair and 
embedding energies. As p tends to 2q the minimum in 



the effective pair potential disappears*^ Therefore, one 
tends to see a progression from potentials that favour 
icosahedra to decahedra to close-packed clusters and fi- 
nally disordered clusters as p — > 2q. This is somewhat 
simplistic since as one goes further away from p = 2q the 
other parameters of the potential play an increasing role 
in the relative depths of the pair and embedding func- 
tions. 

Global optimization of the lead clusters to locate 
the lowest-energy structures was achieved using the 
basin-hopping method^ supplemented by reoptimiza- 
tions of databases of structures found previously for other 
clusters ^4 



III. RESULTS 

Putative global minima up to N = 150 were found. 
For reference, their energies and point groups are given in 
Table[|] The coordinate files for all the global minima can 
be obtained online at the Cambridge Cluster Databases- 
Results up to A=55 have previously been reported by Lai 
et alr^ but the current results improve these at ^=43, 
50, 51 and 54. The size-dependence of the energies are 
plotted in Fig. [3 and a selection of clusters are depicted 
in Fig. 

Of the clusters for TV > 13, 70 are decahedral, 29 are 
close-packed (both perfect face-centred-cubic structures 
and those involving a twin plane), 6 have a distorted dec- 
ahedral form recently characterized for zinc and cadmium 
clusters^ (e.g. 7V=16), three are icosahedral and three 
have polytetrahedral structures. The remaining 27 have 
no assigned structure, and are somewhat disordered. As 
is fairly common, these 'disordered' structures appear at 
the smaller sizes in between the magic numbers at N=13, 
38 and 54 that are associated with ordered clusters. Ex- 
amples at ^=23, 30 and 40 are shown in Fig.|3J 

Interestingly the uncentred Mackay icosahedron at 
iV=54 is more stable than the usual 55-atom Mackay 
icosahedron. The strong compressions at the centre of 
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FIG. 2: Energies of the lead clusters modelled by the (a) Gupta and (b) glue potentials relative to Save, a four-parameter fit to 
the energies. £°" pta = -2. 0631/V + 0.97787V 2/3 - 1.0084iV 1/3 + 1.2917; S|^ c = -2.02517V + 1.6608iV 2/3 + 1.36627V 1/3 - 0.8634. 



the icosahedron resulting from shrinkage of the surface 
bonds makes it unfavourable for this site to be occu- 
pied. This type of vacancy was first noted by Boyer 
and Broughton for very large Lennard- Jones Mackay 
icosahedra, 16 but is not uncommon for much smaller 
metal clusters because of many-body effects. 

For N > 59 all the clusters are either decahedral or 
close-packed. A significant proportion of the close-packed 
clusters are in the region 7V=85-100, where decahedral 
clusters with five atoms along the fivefold axis are becom- 
ing less stable, yet those with six atoms along this axis 
are not yet very competitive. It is noticeable that the 75- 
atom complete Marks decahedron is particularly stable 
(Fig. Pla)). However, the other complete Marks deca- 
hedra with 101 and 146 atoms are not as prominent in 
Flg-GI a ) because of their rather prolate and oblate struc- 
tures, respectively. Indeed, a number of the intermediate 
asymmetric Marks decahedra have similar stabilities. 

The observed structures tally quite well with the value 
of p/q (2.625). There are fewer disordered and close- 
packed structures than gold (e.g. p/q — 2.457 (Ref. 0)), 
but less icosahedral and more decahedral structures than 
for silver (e.g. p/q =3.003 (Ref. 0))i±£ That is, the rela- 



tive weights of the embedding and pair terms is such that 
icosahedra are for the most part not favoured because of 
unfavourable embedding energies associated with com- 
pression of the cluster centre, but decahedra are more 
prevalent than close-packed clusters because of a more 
favourable pair energy due to their greater proportion of 
{111} faces. 

When these results are compared to those previously 
obtained for the glue potential^ there is virtually no over- 
lap. The pattern of magic numbers is completely differ- 
ent (Fig. 0) and none of the structures are the same for 
N > 8. A selection of stable clusters for the glue poten- 
tial are depicted in Fig. 0] The clusters do not exhibit 
any of the usual structural forms, and although most are 
fairly disordered, some strong magic numbers are appar- 
ent in Fig. E^b) • The strongest of these is from the highly 
symmetrical 148-atom hexagonal barrel structure. Also, 
of interest is the 54-atom cluster as it is related to the 
uncentred icosahedron favoured for the Gupta potential 
by a series of twists to the outermost shell around one of 
the fivefold axes. 

The long-range character of p(r) is in part because the 
potential reproduces the small surface energy difference 




FIG. 4: A selection of particularly stable clusters for lead clusters using a glue potential. 
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TABLE I: Energies and point groups of the global minima for lead clusters modelled using a Gupta potential. For N > 13 a 
structural assignment has been made, where i stands for icosahedral, d for decahedral, fee for face-centred cubic, cp for close- 
packed, dd for distorted decahedral and pt for polytetrahedral. If no assignment has been given the cluster has no apparent 
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between the fee {111} and {100} faces of leadi^ Although 
an atom on a {100} faces has one less nearest neighbour, 
this is in part compensated because it has more next- 
neighbours at v2?"o that make a significant contribution 
to p(r). For the clusters this does not lead to {100} 
facets but most of the surfaces of the clusters involve 
some combination of squares and triangles. 

That more disordered structures are favoured for the 
glue potential compared to the Gupta potential is sim- 



ply because of the greater dominance of the embedding 
term. However, the reasons why the glue global minima 
have the best embedding energy is more subtle. In par- 
ticular, it is related to the position of the cutoff, which 
occurs between the second and third neighbour shells of 
the ordered structures. By comparison the disordered 
clusters have a much broader spectrum of pair distances, 
and it is in fact the contributions to the pi values from 
distances near to the cutoff that tip the balance in favour 
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of the observed structures* 2 " The position of this cutoff is 
not physically motivated, but has been chosen for com- 
putational convenience. As the structures do not seem 
to reflect robust features of the potential, they showed 
be viewed with appropriate skepticism. 

Unfortunately, experiments on the structure of lead 
clusters are few. There have been some mass spectro- 
scopic studies of very small clusters, but these are in 
the size range (N < 25)ji2*i*i where empirical potentials 
are expected to be least applicable. The most relevant 
experiments are those by Hyslop et al. using electron 
diffraction. 20 However, for the size range most relevant 
to the present study, no adequate fit to the diffraction 
patterns was possible, whereas the larger clusters appear 
to be dominated by decahedra. 

IV. CONCLUSIONS 

It is the contention of this paper that to make progress 
in predicting the structure of clusters, understanding the 
relationship between the interactions and the resulting 
structure is important. Such an approach will provide a 
better general framework for understanding the criteria 
for a potential to be successful, and specific insights into 
the utility of a particular potential. 

That the two lead potentials we consider suggest a to- 



tally different pattern of structures illustrates some of the 
problems for predicting cluster structure. However, given 
the differences in the potentials when compared in their 
effective pair format, the resulting structural differences 
are not in themselves so surprising. 

The results also illustrate some of the pros and cons of 
using potentials with a specified form as with the Gupta 
potential, versus those with flexible forms. The former 
have less opportunity to produce manifestly unreasonable 
structures, but this is because of their more limited range 
of structural behaviour. If the real structures lie out- 
side of this range, then such potentials have no chance of 
successful structural prediction. For example, the Gupta 
potential predicts that decahedral structures dominate at 
the larger sizes in this study, whereas the experimental 
electron diffraction suggest none of the standard formsiSi 

In contrast, flexible potentials have more chance of pre- 
dicting some of the unusual structures that may occur 
at small sizes. However, there can be unforeseen conse- 
quences in regions of configuration space that lie outside 
of the possibilities considered in the fitting process. In 
the current case, the lead glue clusters are able to ex- 
ploit questionable features of the potential, and so the 
structures are probably not realistic. This is not because 
this potential is a bad one — it has been widely used with 
reasonable success. Rather, it illustrates how predicting 
cluster structures is a very stringent test of a potential. 
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